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Abstract 

We have considered the dynamical evolution of cellular patterns controlled by a stochastic 
Glauber process determined by the deviations of local cell topology from that of a crystalline 
structure. Above a critical temperature evolution is towards a common equilibrium state from 
any initial configuration, but beneath this temperature there is a dynamical phase transition, 
with a start from a quasi-random state leading to non- equilibrium glassy freezing whereas an 
ordered start rests almost unchanged. A temporal persistence function decays exponentially in 
the high temperature equilibrating state but has a characteristic slow non- equilibrium aging- 
like behaviour in the low temperature glassy phase. 

The hexagonal tiling is the best partition of the plane in equal cells. It solves both the 
packing and the covering problems. It is the space-filling assembly of equal cells with the 
minimal interfacial extension. 

It is regular, perfect and beautiful. Surprisingly, however, it is never realized in natural 
biological tissues where a relevant amount of disorder is always present |]1], |2| . This is reminis- 
cent of the situation found in many covalently bonded solids, where, despite the lowest energy 
state being crystalline, in practice amorphous glassy structures are the common quasi-stable 
states ||. Guided by this observation and by recent analysis of idealized infinite-ranged 
spin glasses [Q and in a desire to provide both a linkage to studies of covalent glasses and 
a minimalist model for understanding, we have studied the dynamics of the simplest 'cova- 
lently bonded' network subject to a simple stochastic dynamics characterized by a minimal 
set of control parameters, namely (i) deviations of local topology from crystalline and (ii) a 
temperature. The model has similarities with those of the classic works of Weber et. al.|| 
and Wooten et. al.|| [7], but complements and extends those studies by emphasising (and 
quantifying) not only the freezing of single time measures and slow relaxation, but also the 
apparent non-equilibrating non-stationary (aging) character of two-time correlation observ- 
ables in the glassy phase. It therefore provides a link between the two fields of glass studies 
as epitomized by reference Q on the one hand and references 0, |(J on the other. It should 
provide a valuable starting point for further theoretical analysis of systems complementary 
to those studied in[|J and in recent work on glass transitions in model systems with cen- 
tral forces[Q. It is also extendable, at least in principle, to higher dimensions/more- armed 
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Figure 1: (a) A topologically stable cellular partition (froth), (b) A Tl move. 

vertices. Our study combines Monte Carlo simulation and approximate analysis and con- 
centrates on simple but novel one-time and two-time observables, introduced to probe issues 
of the character of those emphasised in@]. 

In the hexagonal tiling each cell has six neighbours 

and three cells meet at a common vertex. It is a regular two-dimensional, minimally- 
connected pattern. More generally we can consider disordered cellular patterns with 3- 
connected vertices such that the average number of sides per cell is 6 (Euler's theorem ||). 
Locally, the degree of homogeneity can be measured in terms of the number of sides per cell. 
If rii is the number of sides of cell i, then qi = 6 — rii is a measure of its deviation from the 
hexagonal configuration. This quantity is the topological charge of the cell and is related to 
the local curvature. Local vertex re-arrangements conserve the total topological charge but 
allow migration and the annihilation of opposite charges on adjacent cells. When penalties 
are imposed for non-zero charges this mechanism leads the system to self organize. 

Let us now turn to quantification. We characterize any structure by an 'energy' and 
consider random sequential stochastic dynamics driven by changes in that energy. The 
energy is a measure of the inhomogeneity in the structure. In the spirit of minimalism we 
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take it to be given by the 'topological distance' from the hexagonal lattice configuration: 



N N 



E = ^N = Y.{qi) 2 = Y,^-n t f. (1) 



i=i i=i 

We consider a system conserving cell number (N). It evolves by Tl moves, which are 
topological transformations consisting of an exchange of neighbours between 4 cells [see 
fig. (|I|)] . Consider a Tl move performed on a system of four cells where the two adjacent cells 
have n a and n b sides and the two second neighbour cells have n c and n d sides respectively. 
After the Tl move, the cells with n a and n b sides lose one edge each, whereas the cells with 
n c and n d sides each gain an edge. The change of the energy associated with this move is 

AE(n a , n b \ n c , n d ) = 2(2 + n c + n d - n a - n b ) 

We consider dynamics of a Glauber-Kawasaki type, where the probability of this Tl move 
is given by 

Yi(n a ,n b ]n c ,n d ) 



1 + exp(/3AE(n a , n b ; n c , n d )) 
x(l-5 na , 3 )(l-5 nb ,s)(l-6 c>d ) (2) 

This probability allows a dynamical evolution even at zero temperature {(3 = oo) if it reduces 
the energy or leaves it unchanged. In Eq.(|2|) the first two ^-function terms have been intro- 
duced to exclude moves that generate two-sided cells. The last term forbids moves which 
generate cells which are self neighbours, preventing the formation of tadpoles. Thus neither 
two-sided nor single-sided cells can arise if not present at the start. 

Before considering the dynamics, let us consider the predictions of equilibrium. The 
partition function is the sum over all the possible froths with Boltzmann weights. Some 
limitations apply on the accessible values of the rii. 

For instance, the average number of neighbours {(n)) is fixed to be equal to 6 0. Moreover, 
two-sided cells are not admitted (therefore > 2). Self-neighbouring cells are also not 
admitted (which trivially implies ri^ < N), but also other less trivial combinations of ri^ 
might be forbidden. In general, to find the sets of admissible {n{\ which lead to possible 
froths is a very difficult task. 

We approximate the partition function: 

N N-l 

Z(p,\,N)~H £ exp[-/3£(6-n l ) 2 -A5:(6-n i )], 

j=lrij=3 i % 

where A is a Lagrange multiplier fixed by the constraint (n) = 6. Correspondingly, the 
probability of an n-sided cell in the system is 

p(n) ~ p(6) exp[-/3(6 - nf - A(6 - n)]. (3) 

When iV -> oo and (3 = 0, A = ln(3/4) and p(n) = 16/27(3/4) n , with /i 2 = 12. 

To study the dynamics we performed extensive computer simulations on froths with iV = 
100 1 72 (= 3 1 7 2 ) cells and periodic boundary conditions. We started from two different 
systems: (i) a very disordered network, obtained by performing 10 4 iV Tl moves on edges 
chosen at random from an ordered hexagonal seed; (ii) a perfectly ordered pattern (the 



3 



10 
H2 



0.1 





(a) (3 = 0.5 


: A 
A 




: (U)a 




A 





10 



100 t 1000 




1000 



Figure 2: (a) \i<i versus the number of attempted moves in units of N at (3 = 0.5 with starting configurations 
(i) and (ii) (triangles down and up), (b) (3 = 3. The lines are the theoretical predictions for thermodynamical 
equilibrium. 
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hexagonal tiling). In the simulation, 1100 N Tl moves were attempted, with a probability 
given by Eq.(|2]), on edges chosen at random. First we consider a moderate-high temperature. 
In fig. (2a) /i2 versus the number of attempted moves when (3 = 0.5 is shown for the two 
starting froths. These systems evolve from the two initial states (i) and (ii) with /i 2 ~ 13 
and fi2 — [[J respectively, toward a final equilibrium configuration with a common /i 2 — 1 
reached after about lOOA attempted moves. Second we turn to a lower temperature. In 
fig. (2b) /i 2 is plotted versus the number of attempted moves for (3 = 3. In this case, after 
HOOiV attempted moves, the equilibrium configuration is not yet reached, the two patterns 
generated from the two initial states (i) and (ii) are still statistically different, and the 
approach toward a common equilibrium configuration is very slow. 

In fig.(||D p(n) and /z 2 are plotted for several simulations performed at different tempera- 
tures, starting from the two initial states (i) and (ii) and attempting HOOiV moves. At high 
temperatures (/3 _1 > 1) the two final distributions p(n) and /x 2 corresponding to the starting 
configurations (i) and (ii) 

coincide and are in good agreement with the analytical prediction ||. The system has 
reached thermodynamical equilibrium and we designate the state as a 'liquid'. At low tem- 
peratures (P^ 1 < 1) the p{n) and /i 2 reached do not coincide any more. The dynamics is 
slow and the system no longer reaches the thermodynamical equilibrium state in the time- 
scale studied. In the low temperature region for T < 0(1), an ordered start rests almost 
unchanged and a random start seems to lead to a freezing with p(n) close to the form given 
by Eq.(|3p with (3 = 2.4. Hence we identify ~ (2.4) -1 as the dynamical freezing transition 
for the time-scale studied and we designate the corresponding metastable phase as a glass. 

In order to better characterize these two different dynamics at high and low temperatures, 
we have also studied a two time persistence function C(t w + T, t w ) [TT|. It counts the fraction 
of cells that have not been involved in a Tl move between the time t w and t w + r. In 
fig.(^) the persistence functions for j3 = 0.5 and 10 8 , are reported versus r for several value 
of t w . At high temperatures the persistence function appears time-translationally invariant 
and decays exponentially fast as in a liquid, whereas at low temperatures it shows a slow 
aging-like behaviour as in a glass [pi], f|. Thus, at low temperatures the system has not 
equilibrated, even though the single-time measures are no longer evolving significantly. 

Let us now turn to theoretical expectations for C and consider a given cell i. The prob- 
ability that between times t and t + 1 one of its edges or one of the edges incident on its 
vertices (2rij edges in total) is chosen (among the 3N edges in the system) to attempt a 
Tl move is 2rii(t)/3N. This Tl move is effectively performed with a probability Ai(t) that 
depends on the local configuration as in Eq.(^). The probability that a given cell % is not 
involved in any move between the time t w and t w + r is therefore 

*™+ T 2n(t) 

c l (t w + r,t w )= JJ [l - —ILZ^t)] 

t — t w 

~ex P [-£ T ^U(t)] , (4) 

and the expectation value of the persistence is C(t w + r, t w ) = J2i Ci(t w + r, t w ). 

At high temperatures in thermodynamical equilibrium we expect that the local configu- 
ration around the cell i evolves through a set of similar configurations typical of the whole 
system. Therefore, for large r, Yft=t T 2 n i(t)Ai(t) — 2nAr. This leads to a persistence func- 
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Figure 3: (a)^2 versus temperature (T = with starting configurations (i) and (ii) (triangles down 

and up) dynamically equilibrated for 11007V moves. The full curve is the theoretical prediction at the 
thcrmodynamical equilibrium, (b) Side distributions (p(n)) at low temperatures in the "glass phase" (/? = 
3, 10, 100). The curve is the equilibrium prediction for (3 — 2.4. (c) Side distributions at high temperatures 
in the "liquid phase" (/3 = 0,0.5, 1). The full lines are theoretical equilibrium predictions. 
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tion which is independent of t w and decreases exponentially fast. This prediction is in good 

agreement with the results from simulations for (3 < 1 where we obtain best-fit values for 

nA equal to 0.9, 0.83 and 0.42 for [3= 0, 0.5 and 1 respectively. 

At low temperatures, almost only 5, 6 and 7 sided cells are present in the system. We 

can restrict consideration to Tl moves which either leave the energy unchanged or reduce 

it. The former act to move pairs of 7-sided cells, pairs of 5-sided cells and 5-7 couples with 

adjacent pairs of 6-sided cells. 

Pairs of 5-7 couples are brought together by this motion and annihilate to produce four 
6-sided cells. The probability that between time t and t + 1 a 5-7 couple produces a 

Tl move in a given cell i is proportional to the number M 5 | 7 of 5-7 couples in the system: 

nA oc M 5 | 7 . 

The cellular pattern evolves reducing the number of 5-7 couples. This happens when two 
5-7 couples are close together and a Tl move 

annihilates both. The number of couples that annihilate per unit of time ( J 1 ' 7 ) is propor- 
tional to the probability for two couples to be close together which, ignoring correlations, is 
proportional to (^^) 2 - This implies M 5 | 7 oc (t + to)' 1 and consequently J2t=tj ^ n i{t)Ai(t) ~ 
a\n(t w + t + t) — a\n(t w + t ), with a a coefficient which depends on the details of the 
annihilation process. By substituting into Eq.(Q), we get C(t w + r,t w ) ~ Co( t tv ^^ T ) a . In 
fig.® this expression is compared with the results from the simulations for (3 = 10 8 . When 
t w < 10N the behaviour of the persistence function is in good qualitative agreement with 
this crude theoretical prediction with fit parameters a ~ 2.5, to — 2iV and Cq = 1, but for 
large waiting times t w > lOOiV the dynamics become slower. 

In fact the recent spin- glass work[|J emphasises that aging in two-time quantities continues 
for both waiting time t w and relative time r much greater than the characteristic settling 
time-scale for one-time quantities. Hence it would be of interest to extend these studies to 
much longer run-times, as well as to higher dimensions. 

In summary, in a minimalist model of a continuous network subject to stochastic dynamics 
determined by a simple local topological 'energy' and a randomizing temperature we have 
identified a dynamical phase transition and have distinguished two phases: "liquid" at high 
temperatures and "glass" at low temperatures. The dynamics is fast and equilibrating from 
any starting state in the liquid phase and slow and aging in the glass phase which is achieved 
from a random starting state. In our explicit characterization the dynamics become slow for 
(3 < 1 and the transition to the glass phase is around the point j3 = 2.4. 

Although the micro-dynamics we have considered is idealized and chosen for minimalist 
study, it is interesting to compare our results with those observed in nature. Most undiffer- 
entiated biological tissues are space-filling assemblies of cells where the side distribution is 
centered around n = 6 and the second moment [12 takes values between 0.5 and 1.2; for in- 



stance, /i 2 = 0.53 in human epidermis [12], 0.6 < //2 < 1.1 in vegetable leaves fll3 |, ^2 = 0.68 



in cucumber epithelia and [12 = 1.00 in human amnion ]![]. 

We have found slow dynamics in systems with low amount of disorder ((3 > 1 and /i2 < 0.5) 
and fast dynamics in systems with higher disorder (/3 < 1 and fi 2 > 0.5). One might speculate 
that the 'ideal' biological tissue must fit the compromise between low disorder (homogeneity) 
and 

fast dynamics (efficient recovering of perturbations), implying a structure with a value of 
fi2 a little above 0.5. 
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Figure 4: Persistence function C(t w + r,t w ) for t w = (4 )N/3 with k = 1, ..,5. (a) At high temperatures 
(j3 = 0.5) the persistence function is independent of t w and decreases exponentially fast, (b) At low 
temperatures ((3 = 10 8 ) C(t w + r, t w ) depends on t w and shows a slow decay. The curves correspond to fits 

) a with a 



of the form C(t w + r, t w ) ~ C (^^^ r ; 



2.5, t = 2iV and C = AT. 
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